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STABILITY OF CASTERING WHEELS FOR AIRCRAFT LANDING GEARS 

By Arthur Kantrowitz 



SUMMARY 

A theoretical study has been made of the shimmy of 
castering wheels. The theory is based on the discovery of 
a phenomenon called kinematic shimmy. Experimental 
checks, use being made of a model having low-pressure 
tires , are reported and the applicability of the results to 
full scale is discussed. Theoretical methods of estimating 
the spindle viscous damping and the spindle solid friction 
necessary to avoid shimmy are given. A new method of 
avoiding shimmy — lateral freedom — is introduced. 

INTRODUCTION 

In many installations of castering rubber-tired wheels, 
there is a tendency for the wheel to oscillate violently 
about the spindle axis. This phenomenon, popularly 
termed "shimmy," has occurred in some airplane tail 
wheels and has been corrected in two ways: First, by 
the application of friction in the spindles of the tail 
wheels; and, second, by locking the wheels while taxying 
at high speeds. Shimmy is common with the large 
wheels used as nose wheels in tricycle landing gears and, 
because it is impossible to lock the wheels, friction in the 
nose-wheel spindle has been the sole means of correction. 
Because the nose wheel is larger than the conventional 
tail wheel and usually carries a greater load, the larger 
amounts of spindle friction necessary to prevent shimmy 
are objectionable. 

Several experimental investigations undertaken to 
find methods of avoiding shimmy showed that shimmy 
could be prevented in small tail wheels by friction in the 
spindle. The most thorough work on the general 
subject of stability of castering wheels is that of Becker, 
Fromm, and Maruhn (reference 1) on shimmy in auto- 
mobiles. Because they do not include the effect of 
lateral tire deflection, it is believed that their results 
would, in general, not be applicable to airplanes. 

The present paper presents a theoretical and an ex- 
perimental study of the problem of the stability of 
castering wheels for airplane landing gears. On the 
basis of simplified assumptions induced from experi- 
mental observations, a theory of the phenomenon is 
presented. The theory is then compared quantitatively 
with the results of model experiments. 



THEORY OF THE STABILITY OF CASTERING WHEELS 

KINEMATIC (OR STATIC) SHIMMY 

Some preliminary experimental results on shimmy 
were obtained by the N. A. C. A. with the aid of the 
belt-machine apparatus shown in figure 1. This 
machine consists of a continuous fabric belt mounted 
on two rotating drums and driven by a variable-speed 
electric motor. Provision is made for rolling a caster- 
ing wheel up to about 6 inches in diameter on the belt 
in such a way that it is free to move vertically but not 
horizontally. 

On this belt machine, the following phenomenon (see 
fig. 1) was discovered while pushing the belt very 
slowly by hand. With the wheel set at an angle with 
the belt as in (a) and the belt pushed slowly, the 
bottom of the tire would deflect laterally as is shown 
in (b). When the belt was pushed farther, the wheel 
straightened out gradually as is shown in (c). The 
bottom of the tire would then still be deflected, however, 
and the wheel would continue to turn as in (d). The 
wheel would thus finally overshoot, as shown in (e) 
and (f). The process would then be repeated in the 
opposite direction. 

Figure 2 is a photostatic record of the track left by 
the bottom of the tire on a piece of smoked metal. 
Two things will be noticed: First, that the bottom of 
the tire did not skid; and, second, that the places where 
the wheel angle is zero (indicated by zeros on the 
track) correspond roughly to the places where the lateral 
deflection of the tire is a maximum. Thus the wheel 
angle lags the tire deflection by one-quarter cycle. 

It was noticed that the oscillation could be inter- 
rupted at any point in the cycle by interrupting the 
motion of the belt without appreciably altering the 
phenomenon. From this observation it was deduced 
that dynamic forces play no appreciable part in this 
oscillation. 

The distance along the belt required for one cycle 
was also found not to vary much with caster angle or 
caster length (see fig. 3). Caster length was therefore 
considered not to be of fundamental importance in this 
type of oscillation. 

l 
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Figure 2.— Record of the track of a wheel executing a kinematic shimmy. 



It should bo pointed out that, in order to observe 
the kinematic shimmy, lateral restraint of the spindle 
is necessary to prevent the spindle from moving later- 
ally when the bottom of the tire is deflected and thus 
neutralizing the tire deflection. This restraint is sup- 
plied by the dynamic reaction of the airplane when the 
airplane is moving forward rapidly but is not ordi- 
narily present when the airplane is moving forward 
slowly. It has been observed, however, on airplanes 
towed slowly with two towropes so arranged as to 
provide some lateral restraint. 



I II us t rated 
directions 
positive 




^ Spindle axis 



Figure 3.— Definition of terms. 

Figure 1 (a) shows that, when the center line of the 
wheel is at an angle 0 (see fig. 3) with the direction of 
motion, the bottom of the tire deflects. This situation 
is represented schematically in figure 4. It is seen that 
a typical point on the peripheral center line must have 
a component of motion perpendicular to the wheel 
center line if the tire is not to skid. Thus 

d\=^ —sin 6ds 



(The minus sign follows from the conventions used as 
shown in fig. 3.) Since only small oscillations are to 
be considered, the approximation 



d\ 



(1) 



may be substituted. 

The effect of tire deflection on d will now be con- 
sidered. For the purposes of rough calculation, it w ill 
be assumed that, as illustrated in figure 5, the projec- 
tion of the peripheral center line on the ground is a 
circular arc intersecting the wheel central plane at the 
extremities of the projection of the tire diameter. (See 
fig. 1 (d).) Thus, in figure 5, r is the tire radius. (It 
will be assumed for the time being that the caster length 



Direction of motion 



Km. irk 



^Wheel central plane 



„ Path of typical point 
t'l on tire relative to 
if' I wheel = ds 

& I 



Deformed peripheral 
center line 



-Geometrical relationships involved in kinematic shimmy, illustrating 
that d\ = -dds. 



Lateral displace 
ment of typical 
point = dX 





and the caster angle are zero.) Now if the tire is de- 
flected in the form of a circular arc, then the condition 
that the torque about the spindle axis be zero is that 
the strain be symmetrical about the projection of the 
wheel axle on the ground. Clearly, if the wheel is dis- 
placed, it will be turned about the spindle axis by the 
asymmetrical elastic forces until, if it is allowed time 
to reach equilibrium, the symmetrical strain condition 
is reached. Thus, if the tire is deflected an amount X 
as in figure 5 (a) and if the wheel rolls forward a distance 
ds to the condition shown in figure 5 (b), in order for the 
strain to remain symmetrical the wheel must turn 
about the spindle axis an amount dd. From figure 5, 
Rdd=ds. The value of R may be readily obtained 
from geometry in terms of r and X. Thus 22 2 =r 2 + 
(22— X) 2 , from which, if X 2 «r 2 , it is seen that R=r 2 /2\. 
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Then substituting for R, 

dd 



as r 



(2) 



If the caster length is finite, the strain will not be 
symmetrical about the axle, as was assumed here, but 
will be symmetrical with respect to some line parallel to 
the axle but a certain fixed distance ahead. Hence, the 
essential elements of the geometry are unchanged and 
all the reasoning that led to equation (2) is still valid 
for this case. 

Since the phenomena represented by equations (1) 



Wheel central plane 



In order to take account of tires for which the 
assumption made concerning the projection of the 
peripheral center line is not quantitatively valid, an 
empirical constant K x will be used in place of 2\r 2 in 
equations (2) and (3), thus obtaining 



and 



ds 2 1 



(4) 



(5) 



Center line of dis- 
torted periphery 




Projection of peripheral center 
line not in contact with ground 



The constant Ki can be measured by observing the space 
interval of kinematic shimmy. Where experimental 
values of K x are available, they will be used rather 
than the rough theoretical value 2/r 2 . 

DYNAMIC SHIMMY 

In the foregoing derivation for the oscillation 
called kinematic shimmy, it was assumed that 
the strain of the tire was always symmetrical, 
that is, the wheel was moving so slowly that any 
torque arising from dynamic effects involved in 
the oscillation would be negligible. For this case, 
from equation (4) 



Wheel central plane - 



Center line of dis- 
torted periphery 



-de 
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Ids 
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(b) 

,«* Axle and axis of symmetry 



-Ti- 



de 



Form er axis of symm e try 



Figure 5.— Geometrical relationships involved in kinematic shimmy, illustrating that 

dd=2\jr*=K\\. 



and (2) occur simultaneously, they must be combined 
to get the total effect. Thus 



#6 

ds 2 '~ 



(3) 



This differential equation corresponds to a free simple 
harmonic oscillation occurring every time (he wheel 
moves a distance S=2T/^2/r 2 =irr-yj2. 

Measurements of the space interval S of kinematic 
shimmy have been made for three similar tires of the type 
illustrated in figure 1. These tires all had radii of approx- 
imately 2 inches so that the theoretical interval was 
about 0.74 foot. Their experimental intervals were 
0.65 foot, 0.74 foot, and 0.79 foot. This agreement is 
closer than might have been anticipated in view of the 
roughness of the assumption. It will be seen from 
equation (1) that X is one-fourth cycle out of phase 
with d, which is in agreement with the information 
obtained from figure 2. 



If, now, the wheel is assumed to be moving at 
a velocity such that the effect of the moment of 
inertia about the spindle axis is significant, then 
the strain can no longer be symmetrical and, for 
small asymmetries, the torque exerted by the tire 
on the spindle will be proportional to the amount 
of the asymmetry. Thus, the value in paren- 
1 1 L6S6S will no longer be zero but it can be assumed 
that it will be proportional to the dynamic torque ; 
hence 

d 2 9_ ( 1 de\ 
dt'-^V'K.Ts) 



where C x is an appropriate constant of proportionality 
and includes the moment of inertia about the spindle 
axis I w . If the forward velocity V of the wheel is 
constant 

4 2 e 



d 2 6 T T 

dT 2 ~ v 



ds 2 



and 



(6) 



The constant C x may be measured by deflecting 
the bottom of the tire a known amount, moving the 
wheel forward, and balancing the torque T exerted by 
the tire on the spindle so that 0 stays constant. The 
method of deflecting the bottom of the tire a known 
amount will be described later. In this case (6) becomes 



T 
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from which (\ may be found. It has been found that 
Ci increases with increasing caster angle. Thus, for 
a tire like the one in figure 1, d was 71,000 radians 
per second 2 per foot for a caster length L of 0.17 inch 
(caster angle, 5°; no fork offset) and was 104,000 
radians per second 2 per foot for caster length 0.68 inch 
(caster angle, 20°). 

In the study of kinematic shimmy, it was also seen 
that the only change in X was due to the fact that a 
component of the forward motion was perpendicular 
to the central plane of the wheel. This circumstance 
is expressed by equation (1). It is, however, found 
that, if the spindle is clamped at 0 0 and the bottom 
of the tire is deflected, the deflection will gradually 
neutralize itself; that is, the bottom of the tire will 
roll under the wheel. Thus the asymmetrical (dd/ds=0) 



Peripherol > f 

center line \ 



dA-i 



ds 



Contacted 
region 



I 



Wheel 
axle 



Wheel 
central 
plane 



Figuke 6— Illustrating the contribution of asymmetrical strain to dX/ds. 

strain that exists in this case contributes to d\/ds. 
The case of 0 and d6/ds zero is illustrated in figure 6. 
If the effect is again supposed to be proportional to 
the cause, there is obtained 



ds 



= -C 2 \ 



(7) 



In this equation, the constant C 2 is a geometrical 
constant of the tire that can be obtained from static 
measurements. The order of magnitude of C 2 can be 
obtained by assuming that the periphery of the tire 
intersects the extremity of the extended central plane 
of the wheel. In that case C 2 =l/r. 



If 6 is not zero, there will be a component of the for- 
ward motion contributing to d\/ds. As in equation 
(1), this component will be —6. Adding this compo- 
nent to the part of d\/ds due to asymmetry, then 
(dd/ds still assumed zero) 



T~ — — @ — OoX 

ds 



(8) 



de 



This equation expresses that, for g^=0, the contribu- 
tion of the asymmetrical strain to d\/ds was — C 2 \. Also 
for the symmetrical strain, in which case \~ ^J^ S = 0 

(kinematic shimmy), there was, of course, no contribu- 
tion due to asymmetry. Assume now a linear interpo- 
lation between these two limiting cases. Thus, finally, 



ds 



( x ~i;§) 



(9) 



A method of measuring C 2 is provided by equation 
(7) which, when integrated, gives 

If, with the spindle clamped at 0 = 0, the tire is 
deflected a known amount X 0 and rolled ahead a known 
distance and the new X measured, C 2 may be computed. 
It was seen earlier than C 2 was of the order of magni- 
tude of 1/r; that is, it would be of the order of 6 for a 
2-inch-radius tire. The constant C 2 was measured on 
two model tires under different loads and found to be 
G.2 and 3.4. Considerable variations of this constant 
with tire pressure and load have been found. 

A method of obtaining the constant known X neces- 
sary for the measurement of C x is provided through 
equation (8). Here it is seen that, if the wheel is 
pushed along at a constant angle 0, X will increase 

(negatively) until d= — C 2 X, in which case j^=0 and 

equilibrium is reached. 

When the wheel is moving ahead at a finite constant 
velocity, the phenomena represented by equations (6) 
and (9) occur simultaneously. Therefore, to get the 
total effect, combine the two equations, thus obtaining 

<.o) 

The solution for the natural modes of motion repre- 
sented by equation (10) is 

0=Ae a i'+Be a *+Ce a *' C 11 ) 

where the as are the three solutions of the so-called 
auxiliary equation 

c?« 3+ (i; + ? 1/2 > 2+1=0 (12) 
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One of these a's is real and negative and corresponds 
to a non oscillatory convergence. The other roots are 
conjugate complex numbers and correspond to the 
shimmy under consideration. The roots will be of the 
form (/ dzcoi. If the divergence a is positive, the oscilla- 
tion will steadily increase in amplitude (while its ampli- 
tude is not large enough for skidding to occur); and, 
if a is negative, the oscillation will steadily decrease in 
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Figure 7.— Theoretical values of divergence, frequency, and phase angle against 

velocity. 

amplitude and eventually disappear. The meaning of 
the quantity "divergence" may be illustrated by saying 
that it is approximately equal to the natural logarithm 
of the ratio of successive maximum amplitudes to the 
distance between them. The quantity co is equal to 2t 
times the number of oscillations per foot. The fre- 
quency therefore is a>V/2w. The phase angle is obtained 



by substituting for 6 in equation (6) the value obtained 
from the foregoing procedure and solving for X. 

The divergence, the frequency, and the phase rela- 
tions thus derived for typical model tire constants are 
plotted in figure 7. For small velocities (0 to 6 feet per 
second) the frequency corresponds to kinematic shimmy ; 
it is proportional to velocity. The divergence increases 
rapidly, however, because the spindle angle lags on 
account of the moment of inertia about the spindle axis, 
thus allowing more lateral deflection than would occur 
in a kinematic shimmy. On the next half cycle, a 
larger spindle angle is reached and the process repeats. 
As the velocity is further increased, the lag, and hence 
the asymmetry of the strain, further increase until the 
strain becomes almost entirely asymmetrical. For this 



condition, ^ ^-<<X. 
as ill 



Then the restoring torque on 



the spindle is approximately proportional to X. (See 
equation (6).) The tire deflection X (measured nega- 
tively) will, however, still lag somewhat behind 6 be- 
cause, after the wheel is turned through a given angle, a 
certain forward distance is required for equilibrium tire 
deflection to be reached. Thus, the restoring force will 
again lag the displacement. As the velocity increases 
in the high-velocity range, the frequency stays nearly 
constant (see fig. 7) and the distance corresponding to 
a single oscillation increases. Hence this constant lag 
becomes a smaller part of the cycle and the divergence 
decreases at high velocities. 

It will be appreciated that the foregoing theory 
considers only the fundamental phenomena taking place 
in shimmy. Other phenomena occurring simultane- 
ously have been neglected. Some of the more impor- 
tant of the neglected phenomena are: 

1. Miscellaneous strains (other than lateral tire 
deflection) occurring in the tire. A rubber t ire being an 
elastic body will distort in many complicated ways 
while shimmying. In particular, there will be a twist 
in the tire due to the transmission of torque from the, 
ground to the wheel. 

2. Two effects will cause the stiffness constants ol 
the tire to change with speed. First, centrifugal force 
on the rubber will make the tire effectively stiffer at 
high speeds. Second, much of the energy used to 
deflect the tire will go into compressing the air. The 
compressibility of the air will change with the speed of 
compression owing to the different amounts of heat 
being transferred from it. 

3. There will be a gyrostatic torque about the spindle 
axis caused by the interaction of the rotation of the 
wheel on the axle and the effective rotation of part of 
the tire about a longitudinal axis on account of the 
lateral tire deflection. This torque will later be shown 
to have a noticeable effect on the results. 
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(a) Caster angle, 5°. (b) Caster angle, 20°. 

FIGURE 8. — Experimental values of divergence and frequency of shimmy against velocity as determined on belt machine. 



The inclusion of items 1 and 2 in the theory would 
obviously be very difficult. It is therefore necessary 
bo resort to experiment to determine whether the present 
theory gives an adequate descript ion of the phenomena. 
If so, the omission of these and any other items will be 
justified. 

An experimental check on the theory w T as obtained 
by measuring the divergence and the frequency of the 
shimmy on the bell machine at two caster angles and 
at a series of velocities. These measurements were 
made by placing a lighted flashlight bulb on a 6-inch 
sting ahead of a model castering wheel with a ball- 
bearing spindle and then taking high-speed moving 
pictures of the flashlight bulb with the wheel free. The 
photographs were made with time recordings on the 
film, and the belt carried an object that interrupted 
light from a fixed flashlight bulb and thus recorded the 
belt speed on the film. 

The divergence and the frequency of the shimmy 
were obtained by measuring the displacements and 
the times corresponding to successive maximum ampli- 
tudes (while the amplitude was still small enough to 
make all the assumptions valid). The results are 
plotted in figure 8. 

In order to compare these results with the theory, 
the constants C\ } C 2 , and Ki were measured on the same 
tire at the two caster angles by the previously described 
methods. It was found that (? 2 =6.2 feet _1 ; that K x = 
62.5 feet " 2 ; and that, for 5° caster angle, C x = l 1,1 00 
feet "'-second ~ 2 and, for 20° caster angle, C, = 104,000 
feet "^-second ~ 2 . The roots of equation (12) were 
then found and (lie divergence and the frequency of tlx 1 

205891—40 2 



shimmy were computed for a series of velocities and at 
caster angles of 5° and 20°. These results are plotted 
in figure 9 and, for purposes of comparison, the experi- 
mental curves are also reproduced. 

20r 




\ o 
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16 24 
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Figure 9.— Divergence and frequency of shimmy on belt machine compared with 

theory. 
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The agreement between theory and experiment is 
considered satisfactory as regards qualitative results. 
It will be noticed, however, that the theoretical values 
of the divergence are decidedly too large at high veloci- 
ties, say 25 feet per second. A velocity of 25 feet per 
second for the model would, however, correspond to 100 
feet per second for even a small airplane, such as the 
Weick W-1A (reference 2). It was thought that the 
discrepancy might be attributed to a gyroscopic torque 
(mentioned earlier) that was due to the interaction of 
the rotation of the wheel on its axle with lateral motion 
of the bottom of the tire. This lateral motion may be 
considered to be equivalent to a rotation of part of the 
tire about an axis perpendicular to the spindle and the 



20 




\ O 8 16 34 32 40 48 
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Figure 10.— Divergence and frequency of shimmy on belt machine compared with 
theory including gyrostatic torque. 



wheel axle. This interaction gives rise to a spindle 
torque, which must be included in equation (6). The 
numerical value of this torque will be I a ^<p f where l a is 
the moment of inertia about the wheel axle of the part 
of the tire that engages in the lateral motion (say, 
possibly one-fourth of the total moment of inertia of the 
tire); ^ is the angular rotation about the longitudinal 
axis (\f/ = \/r)] and <p is the angular rotation about the 
axle (<p = V/r). Thus, the value of the gyrostatic torque 

is Ia\V/r 2 ; but *=V^ so > finally, the torque is 
I Y-.d£m The angular acceleration due to this torque is 

1* r/Q 



^T^rfX Thug (6) becomeg 

^ 2 d 2 d I a V 2 d\_ n ( x dd/ds\ 
V ds 2 I w r 2 ds~ U \ K K x ) 
Combining this equation with (9) 

v 2 d?e /c H72 c 2 v 2 i a i ye h v 2 de R _ 

Here the constants d, C 2 , and K x are the same as 
before. For the assembly on which the measurements 
plotted in figure 8 were made, I w was 1.06 X 10~ 4 slug- 
feet 2 , I a was 1.09 X 10" 5 slug-feet 2 , and r was 0.165 foot. 
These values were used to compute new theoretical 
divergences and frequencies, which are compared with 
the experiments in figure l 0. It will be seen that, when 
the gyrostatic torques are included, the agreement of 
the theoretical divergences with experimental values is 
considerably improved at high velocities. 

The gyrostatic torque is proportional to d\/ds. 
Hence, at high velocities, it is nearly proportional to 
dS/ds since —X and d are nearly in phase. Thus, this 
torque is approximately in the right phase to damp the 
oscillation. 

EFFECT OF MOTION OF AIRPLANE 

When a shimmying wheel is attached to an airplane, 
the airplane will respond to a certain extent to the 
forces supplied by the shimmy. This lateral motion of 
the nose of the airplane will in turn affect the shimmy 
and does, in fact, reduce the divergence at low veloci( ies. 

In order to calculate an approximation to this effect, 
the airplane will be assumed to be a rigid body. By 
this assumption, any deflection occurring in the struc- 
ture or at the rear tires is neglected. The gyrostatic 
torque due to lateral tire deflection is also omitted 

The following symbols are used in the calculations: 

I z moment of inertia about vertical axis through 

center of gravity of the airplane. 
m mass of airplane. 

E\ lateral force exerted by nose-wheel tire per 
foot deflection. 

6 angle between central plane of wheel and aver- 
age direction of motion. 

\p angle between longitudinal axis of airplane 
and average direction of motion. Measured 
positive in the same direction as 6. 

y lateral distance between center of gravity of 
airplane and path of average mot ion. Posi- 
tive in same direction as X. 

Y lateral force exerted by the rear wheels on the 
airplane. 

l 2} U horizontal distance from nose wheel and pear 
wheels, respectively, to center of gravity. 
l=h+l 2 . 



The side force developed by the shimmy is E\\ 
The lateral acceleration is 
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Equation (6) is unaffected. Combining these equations 



~ y ~ v ds 2 



in 



(13) 



Similarly, the angular acceleration of the airplane is 
given by 

-i 2 ffx\+u- = ^ =1/2 <g (14) 

There is, in addition, the geometrical relationship 

I— *-i^-*>i t») 

The lateral motion of the nose per unit forward motion, 
assuming that the rear wheels do not skid, is 

This motion will naturally affect the rate of deflection 
of the nose wheel d\/ds; equation (9) becomes 

§=-'-<^e$)+*+^ e» 




(17) 



This differential equation has two real roots, which 
are of little interest to the problem at hand, and two 
complex roots, which give the divergence and the fre- 
quency of the shimmy under consideration. 

An experimental check on the theory was obtained by 
running on a flat surface (see fig. 11) two models similar 
to the one shown in figure 12. The model was launched 
by a compressed-air catapult, which gave it a very 
smooth start. The nose wheel was locked during the 



Camera — 



,Shuiter 



Synchronous mo for 



Figure 12.— Model used for free-model test. 



Compressed air 



Cafapulf 



Flashlight lamp 
Batteries^ / ^Ball bearings 





Valve 



TT 



TJ 



Figure 11.— Apparatus for conducting free-model experiments (elevation view.) 
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acceleration so that the shimmy would not start to 
build up before constant velocity was reached. The 
wheel was freed by a lever tripped by a wire. It then 
came into (lie field of view of an elevated camera and 
the flashlight bulb made a track of light on the camera 
film. The light entering the camera was interrupted 



cases are small. The theory therefore seems adequate 4 
to account for the shimmy of freely rolling models when 
corrections are made for the response of the model to 
the shimmy. 

Full-scale experiments similar to the ones with the 
model would be valuable in fur! her checking the the< »ry. 




Figure 13.— Record from free-model apparatus. 



60 times a second by a shutter driven by a synchronous 
motor. A typical trace is shown in figure 13. The 
curve was read only in the region between the two 
arrows. Several such traces were made at various 
velocities and the divergence and the frequency of the 
shimmy were measured and plotted against velocity in 



One very rough check was made on the W-lA airplane 
( reference 2). The friction in the spindle of this air- 
plane was considerable but, when the friction was re- 
duced, the wheel shimmied. Motion pictures were 
taken of the shimmying wheel and some of the enlarged 
frames are shown as figure 15. The shimmy diverged 
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(a) Model 1. (b) Model 2. 

Figure 14. — Divergence and frequency of shimmy of nose-wheel models comparcdwith theory. 



figure 14. Two models with different mass charac- 
teristics and different caster angles and tires were used 
to get two checks on the theory. The theoretical 
divergences and frequencies in these models are also 
shown in figure 14. 

It will be noticed that both the theoretical and the 
experimental divergences are somewhat smaller at low 
velocities than in the case of the belt-machine shimmy. 
In general, however, the differences between the two 



and finally reached an equilibrium amplitude. It will 
be seen that deflection of the bottom of the tire plays a 
prominent part in the shimmy. The frequency of the 
shimmy of this airplane was measured at a velocity of 
47 feet per second and found to be 9 cycles per second. 
The frequency was then computed for this velocity 
(on the basis of constants found by the method to fol- 
low) and a frequency of 10.5 cycles per second was 
obtained. These results are taken to be indications 
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that the theory will be at least approximately valid for 
full scale. 

In the application of the theory, it is necessary to 
know the constants of the tire to be used. In the ab- 



length; n for the model tires was 0.69 at a caster angle 
of 5° and 1.00 at a caster angle of 20°. 

These expressions for C u C 2 , and K 1 will be used for 
numerical calculations in the rest of the report. 




Figure 15.— Shimmy in W-lA airplane, 64 frames per second. 



sence of such measurements for full-scale tires, it will 
be necessary to use the theoretical values of the param- 
eters K x =2lr 2 and <7 2 = l/r. For I w C l (the spindle 
torque per unit lateral tire deflection for a completely 
asymmetrical strain), the expression nrE x will be used, 
where E\ is the lateral force per unit deflection, r is 
the radius of the tire, and n is an empirical nondimen- 
sional factor to be oh (.-lined from experiments. The 
factor n will, of course, increase with increasing caster 



METHODS OF AVOIDING SHIMMY 

SPINDLE DAMPING 

In the past the only method that has been used to 
avoid shimmy in single swiveling wheels has consisted 
in the application of friction or hydraulic damping to 
the spindle. It will therefore be desirable to set up a 
criterion for the amount of spindle damping necessary to 
avoid shimmy. 
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Iii this calculation and also in later calculations of 
oilier methods of avoiding shimmy, the damping effect 
of gyrostatic torque due to lateral tire deflections will 
be neglected as well as that due to airplane motions; 
these omissions result in small errors and usually give 
a conservative result. 

First, the effect of idealized viscous damping applied 
to the spindle will be calculated, that is, the effect of 
a torque proportional to the angular velocity of the 

spindle, say K® = The angular acceleration 

resulting from this torque KV^/I W must be added to 
equation (6). 



V ds 2 



KVdd 
' I w ds 



(18) 



Combining this equation with equation (9), which is 
unaffected by the addition of spindle friction, 



(7, <fe 3+ 



(vk ± a \ 



In the case of a self-excited oscillation represented by 
a third-order differential equation, Month's discriminant 
provides a simple method of deciding whether the 
oscillation will or will not diverge. If all the coefficients 
are positive (as they are here), then the oscillation will 
diverge if and only if the product of the coefficients of 
the third- and zero-order terms is greater than that of 
the first- and second-order terms. Hence for the shim- 
my to have zero divergence 

or, rearranging and solving for K/I w , 



.VK, 



(20) 



for the condition of zero divergence. 

Only the positive value of .PC is of interest because the 
application of negative damping to the spindle would 
make the wheel unstable. From equation (20) a maxi- 
mum positive value of If is found at the velocity -y/CjC^K h 
which for the model wheels at 5° caster is 13.6 feet per 
second. Thus, the amount of spindle damping neces- 
sary to avoid shimmy is seen to be a maximum near the 
velocity at which maximum divergence occurs with an 
undamped wheel. (See fig. 9.) This velocity will be 
denoted by F 0 . The damping required at this velocity, 
i. e., the maximum damping required, is 

K maI =I,J^( \l 'l+gj-l) a- ^ (2D 

Y K Y \ V 6 2 / radians per sec v 

Equation (21) shows that it is desirable to make the 
caster length small for in this way C u and hence K max , 



may be made small. A rough approximation to K max 
may be obtained from equation (21) by the use of the 
theoretical formulas for the tire parameters given 
earlier. Using these formulas and the value of n 
obtained for 5° caster angle, 



K maT =0A3rjl w E x 



Ib-ft 



radians per sec 

Using these formulas, also 

V 0 =0.59r 2 ^EJI^ c fps 

for the velocity at which this maximum damping is 
required. For the W-lA airplane, using r=0.67 ft, 
7^=0.11 slug-ft 2 (estimated), and E x = 2,400 lb/ft 
(from rough measurements). 

K max =4:.6 lb-ft per radian per sec, and 
y 0 =39 fps 

For the Hammond Model Y airplane, which has had 
considerable difficulty with shimmy, r=0.S3 ft, I w 
(estimated) = 0.33 slug-ft 2 , E x (estimated)=4,130 lb/ft, 
w hich gave 

K max = 13.2 lb-ft per radian per sec 
y 0 =46 fps 

(If full-scale tire-constant measurements as yet 
unpublished are used, equation (21) becomes 

Kmax= constant XrJT w — r lh ~ ft 

v w radians per see 

where the constant may vary between L6 and 32 de- 
pending on the type of tire.) 

The amount of damping required increases with the 
fourth power of the size of the airplane, as would be 
expected. It is necessary in the application of a prac- 
tical viscous damper to insure that adequate damping 
is obtained at all amplitudes; that is, it would be 
necessary to allow for the nonlinearity of the damper. 

Most practical shimmy dampers have used solid 
rather than viscous friction. Solid friction, however, 
is difficult to deal with analytically in a straightforward 
manner. Therefore, resort will be had to the following 
subterfuge. It will be assumed that, if solid friction 
lakes an amount of energy out of an oscillation equal 
to the energy taken out by a given amount of viscous 
damping, it is equivalent to viscous damping for pur- 
poses of damping the oscillation. The energy taken out 
by viscous damping in one-quarter cycle is 

f° ma 'Kd de= r m "KV^d$= \KVa>e max * 

where 0=0 maz sina>.9. The energy taken out by a solid 
friction torque T is T Q max . Therefore, to take out the 
same amount of energy 



T 



(22) 
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(Velocity enters because the independent variable is 
distance.) 

Equation (22) gives an interesting explanation of the 
fact that shimmy (and other self-excited vibrations) 
will usually not start in the presence of solid friction 
until an initial displacement is supplied. From the 
theory of small oscillations the viscous friction required 
to damp a self-excited vibration is independent of the 
amplitude. (See equation (21).) It will be seen from 
(22) that the solid friction required to damp a self- 
excited vibration will be directly proportional to the 
initial amplitude of the oscillation. If solid friction is 
used to damp shimmy, it will be desirable to insure that 
the wheel is centered when it first makes contact with 
the ground. This precaution has been found useful on 
airplanes. 

The value of co 2 for the case of equilibrium viscous 
damping (zero divergence) is equal to the ratio of the 
coefficient of the dd/ds term to the d*d/ds z term in 
equation (19). This fact may be readily verified by 
substituting a simple harmonic solution for 0 and equat- 
ing the sum of the imaginary terms to zero. Thus 



(23) 



It is difficult to obtain a good estimate of the maxi- 
mum angle that should be used in computing the 
amount of solid friction necessary to damp shimmy on 
an actual airplane because the value depends on how 
much energy will be supplied to the shimmy after the 
tire starts to skid. If it be supposed that the energy 
supplied to the shimmy at angles greater than that at 
which skidding starts is entirely taken out by the solid 
friction, then the value of 6 for which skidding begins 
can be put into (22) and a value of the friction will be 
obtained such that shimmy with all initial angles will 
be damped out. This supposition is probably true for 
small caster lengths. The tire will start to skid when 
the lateral force is large enough to overcome the fric- 
tion, that is, when E\\=^fiW 9 where \i is the coefficient 
of friction and W is the load on the swiveling wheel. 
Now, if the wheel moves in a direction at an angle 6 
with the central plane of the wheel, the equilibrium 
deflection of the tire will be given by 

from equation (8). In this way a rough estimate can 
be obtained of the angle 0 O at which skidding starts. 
Thus 

1*1"^ (24) 

Substituting in (22) from (2 4) for 0 max , the value from 
(23) for co, V 0 or -y/CjC^Ki for V, and the value of 
*: mai from (21) for K, then 



(V^i- 1 )" (25) 



is obtained as the amount of solid spindle friction re- 
quired to damp shimmy at the velocity requiring most 
viscous spindle damping (which velocity will probably 
require the most solid spindle friction). The validity 
of formula (25) is restricted to small caster lengths; 
this restriction should not hamper its application be- 
cause the caster length is preferably small. Using the 
theoretical formulas for the tire parameters, there is 
obtained the simple formula (again assuming small 
caster length), 



T max = 0. 17 txrW lb-ft 



(26) 



It should be remembered that the theoretical for- 
mulas for the tire parameters used in the derivation of 
(26) are approximately valid only for tires similar to 
the model tires used in the experiments. Thus (26) 
and all other formulas based on the theoretical tire 
parameters will be approximately valid only for air 
wheels or low-pressure tires. For high-pressure tires, 
it seems likely that the effective r would be smaller 
than the geometrical radius of the tire. It should be 
remembered that no conservative assumptions were 
made in the derivation of (26). 

If now />t=0.55, a value of 19 foot-pounds spindle 
friction is required for the PP-lA airplane with W=S10 
pounds (including propeller thrust) and 44 foot-pounds 
for the Hammond Y airplane with W=n7Q pounds. 
The amount of spindle friction actually present in the 
spindle and steering mechanism of the W-lA was 
measured and found to be of the same order of magni- 
tude as the computed value. This result explains why 
no shimmy was found with that airplane. As men- 
tioned previously, when this solid friction was reduced, 
shimmy occurred in the castering front wheel of this 
airplane. 

(Later full-scale tire-constant measurements as yet 
unpublished indicate that a value T max =0.06 firW lb-ft 
would be more accurate.) 

It will be seen that the amount of solid spindle fric- 
tion necessary to damp shimmy in nose wheels is con- 
siderable even at small caster lengths. As is the case 
for the hydraulic damper, it increases with the fourth 
power of the size of the airplane. 

LATERAL FREEDOM 

The divergence of the shimmy has been seen to 
depend on the interaction of the lateral tire deflection 
with the angular rotation of the spindle. If the cou- 
pling between these two degrees of freedom could be 
sufficiently reduced, the shimmy would not diverge. 

One method of reducing the coupling is to introduce 
a new degree of freedom into the oscillation. Thus, if 
the wheel is free to move in the direction of the axle, 
the tire deflections will depend on the inertia reactions 
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of the wheel and of any other members moving later- 
ally. The maximum mass that the wheel could have 
and still reduce the coupling sufficiently to avoid 
shimmy is calculated in the appendix. The masses of 
wheels used in practice an* found to be small enough 
to prevent shimmy if lateral freedom were used. The 
amount of lateral freedom required is of the order of 
the distance that the wheel will deflect laterally with- 
out skidding. (See the appendix.) 

This system was tried on the belt machine and on the 
W-lA airplane. Shimmy could be prevented if the 
spindle were prevented from turning through an angle 
large enough to make the tire skid. Experiments 
made here and elsewhere without this angular restric- 
tion have failed to prevent shimmy. The necessity 
of limiting the angular travel of the spindle indicates 
that another type of oscillation occurs at large 4 spindle 
angles. It seems likely that this oscillation is due to 
the decrease of skidding friction with increase of rela- 
tive skidding velocity. As is well known, this effect 
is apt to cause a chattering type of oscillation and, in 
the present tests, probably fed energy into the shimmy. 

CONCLUDING REMARKS 

The phenomenon called kinematic shimmy seems to 
he the fundamental phenomenon in the oscillation of 



model wheels. A dynamical theory based on this dis- 
covery gives an adequate description of the shimmy of 
small model wheels. It seems unlikely that shimmy 
of full-scale airplanes will be fundamentally different. 

On the basis of this theory, a method has been pro- 
vided of estimating numerically the spindle friction 
necessary to damp shimmy. An objection to this 
method of avoiding shimmy lies in its indicated inter- 
ference with the -leering of nose wheels. 

A new method of avoiding shimmy by the use of 
la t (Mai freedom for the wheel is introduced but requires 
further investigation before complete design data can 
be given. 



Langley Memorial Aeronautical Laboratory, 
National Advisory Committee for Aeronautics, 
Langlet Field, Va., August 11, 1937. 
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APPENDIX 



LATERAL FREEDOM 

If the lateral acceleration of the wheel is x and the 
mass of the wheel and tire (and other members moving 
laterally) is m', then 

or 

JE\\ dx . 

This lateral motion of the wheel must he added to 
equation (9), 

d\ n(\ dd/ds \ dx 

or 

Tlie effect on equation (6) will be small because only 
small lateral motions will be involved. Combining 
equation (6) with equation (28), 

V 2 <P6 ( q lWfl. gx d6 
Ct ds s± \ v CVKj ds f+ C I m' ds 

It is seen that a damping term has been added to the 
final equation of motion. (Cf. equation (10).) If this 
term is sufficiently large, it will counteract the diver- 
gence of the oscillation. Applying Routh's discrimi- 
nant to obtain the condition of zero divergence and 
simplifying, 

, C 2 E X ■ 

It is noted that m 0 ' is independent of velocity. If m' 
is smaller than m {) ' the shimmy will converge; if greater, 
it will diverge. Putting in the theoretical formulas for 
the constants 

m 0 / =1.45^slugs=47p 2 pounds 

For the W-lA airplane this value would be 11 pounds 
and for the Hammond Y, 22 pounds. Fortunately, 
these values are of the same order of magnitude as the 
weights of the wheels used on these airplanes. 

The necessary amount of lateral freedom will now be 
estimated. The distance is of the order of magnitude 
of the amplitude of the oscillation in r, which corre- 
sponds to a violent enough oscillation to skid the tire. 
After skidding starts, lateral freedom would not be 



expected to have any beneficial effect. For skidding 
to start, 

But 

m'x=--E\\ 

Hence 

Suppose the value of m' used is that for equilibrium 
damping (m 0 ' = C 2 Ex/Ci), then the maximum possible 
xis 

fiWOt 
XmaT ~ExC 2 

But 

*^maar == V W X max 

Hence 

Xmax -V 2 a> 2 ExC 2 

As before, for equilibrium damping, co 2 equals the ratio 
of the coefficient of the dd/ds term to that of the d^Qjds* 
term in equation (29), for m' = C 2 Ex/Ci, 

9 Oj 

and 

%max — 

It will therefore be seen that the minimum amount of 
lateral freedom necessary on either side of the wheel is 
equal to the amount that the tire will deflect laterally 
before skidding when the weight of the wheel is the 
required value. This approximation is rough and no 
conservative assumptions were made other than the 
one implicit in the neglect of whatever spindle friction 
may be present. For the W-l A airplane, the calculated 
amount of lateral freedom required is 0.85 inch on 
either side of the wheel and, for the Hammond Y air- 
plane, 0.91 inch. 

In the application of lateral freedom to an airplane 
wheel, means for keeping the wheel centered in the 
lateral travel must be provided. It can be shown that 
the application of spring-restoring forces sufficiently 
large to accomplish this purpose would tend to make 
the wheel shimmy. It was found, however, in a lateral- 
freedom arrangement that, when the wheel was tilted, 
it would tend to ride on the lowest portion of the axle. 
From figure 16 the point A has a smaller rolling radius 
than the point B but the point A must make as many 
revolutions as point B. Hence, if no skidding occurs, 
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point A will tend to lag behind point B. Then the re- 
sultant couple will turn the wheel about the spindle 
axis in such a direction that the forward motion will 
supply a side force causing the wheel to slide "down- 
hill" on the axle. This effect is found to be more than 



<L wheel 




Figure 16.— Tilted axle. 

large enough to overcome the component of the weight 
of the airplane that would tend to make the wheel slide 
"uphill." If an axle that is convex downward is used, 
the wheel will therefore always seek the lowest point, 
which in this case will be the center of the lateral 



travel. Both theoretically and experimentally, the 
lowest point has been found to be a point of stable 
equilibrium. 

Two methods of applying lateral freedom to a caster- 
ing wheel were tested on the belt machine. The first 
method was simply to allow the wheel to slide along the 
axle. The second method was to allow the wheel to 
rotate about two axes in a vertical plane perpendicular 
to the axle. Shimmying in models could be avoided by 
either method. The second method, however, did not 
involve making the fork any wider. 

Additional full-scale tests of the first system were 
made of the W-lA airplane. An axle curved down- 
ward was used to keep the wheel centered. A radius 
of curvature equal to 6 times the tire radius was arbi- 
trarily chosen. The wheel centered satisfactorily and, 
on a rough turf-covered flying field, shimmy was pre- 
vented by the use of one-half inch of lateral freedom on 
either side of the wheel. On concrete with a lateral 
freedom of three-fourths inch, the wheel could be made 
to shimmy by taxying over a severe bump at a moder- 
ate speed and at an angle that would introduce a large 
initial angular deflection of the wheel. In order to 
prevent shimmy under these conditions, the rotation of 
the spindle was limited by stops to ±13°. The pilot 
did not regard this restriction of spindle freedom as 
objectionable for the ground handling of the airplane. 
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Angle of set of control surface (relative to neutral 
position), 8. (Indicate surface by proper subscript.) 



D, 

V, 

P/D, 

V, 

% 

T, 

Q, 



Diameter 
Geometric pitch 
Pitch ratio 
Inflow velocity 
Slipstream velocity 

Thrust, absolute coefficient C T = 



4. PROPELLER SYMBOLS 



Torque, absolute coefficient C 0 = 



pnW 1 
. Q 



C s , 

v, 
n, 



Power, absolute coefficient C P = pri z]y> 

Speed-power coefficient =*W-p^2 
Efiiciency 

Revolutions per second, r.p.s. 

Effective helix angle = tan -1 ^ 0 \ 

\2Trrny 



~ P n 2 ]y> 

5. NUMERICAL RELATIONS 



1 hp. = 76.04 kg-m/s=550 ft-lb./sec. 
1 metric horsepower— 1.0132 hp. 
1 m.p.h. = 0.4470 m.p.s. 
1 m.p.s. = 2.2369 m.p.h. 



1 lb. -0.4536 kg. 

1 kg=2.2046 lb. 

1 mi. = 1,609.35 m=5,280 ft. 

1 m=3.2808 ft. 
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